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Introduction 

One major obstacle to human space exploration is the possible limitations imposed by the 
adverse effects of long-term exposure to the space environment. Even before human spaceflight 
began, the potentially brief exposure of astronauts to the very intense random solar particle 
events (SPE) were of great concern. A new challenge appears in deep space exploration from 
exposure to the low-intensity heavy-ion flux of the galactic cosmic rays (GCR) since the 
missions are of long duration and the accumulated GCR exposures can be high. Because cancer 
induction rates increase behind low to rather large thicknesses of aluminum shielding, according 
to available biological data on mammalian exposures to GCR like ions, the shield requirements 
for a Mars mission are prohibitively expensive in terms of mission launch costs. Therefore, a 
critical issue in the Human Exploration and Development of Space enterprise is cost effective 
mitigation of risk associated with ionizing radiation exposure. 

In order to estimate astronaut risk to GCR exposure and associated cancer risks and 
health hazards, it is necessary to do shield material studies. To determine an optimum radiation 
shield material it is necessary to understand nuclear interaction processes such as fragmentation 
and secondary particle production which is a function of energy dependent cross sections. This 
requires knowledge of material transmission characteristics either through laboratory testing or 
improved theoretical modeling. Here ion beam transport theory is of importance in that testing of 
materials in the laboratory environment generated by particle accelerators is a necessary step in 
materials development and evaluation for space use [1]. The approximations used in solving the 
Boltzmann transport equation for the space setting are often not sufficient for laboratory work 
and those issues are a major emphasis of the present work. 

In space radiation transport, the energy lost through atomic collisions is treated as 
averaged processes over the many events which occur over even relatively small dimensions of 
most materials and is referred to as the continuous slowing down approximation. It is reasoned 
that the few percent energy fluctuation in energy loss has little meaning for ions of broad energy 
spectra and especially in comparison to the many nuclear events for which uncertainties are still 
relatively large. In contrast, the laboratory testing of potential shielding materials uses nearly 
monoenergetic ion beams in which the interpretation of the interaction with shield materials 
requires a detailed description of the interaction process for comparison to detector responses [2]. 
The development of a Green's function approach to ion transport facilitates the modeling of 
laboratory radiation environments and allows for the direct testing of transport approximations of 
material transmission properties. For a number of years this approach has played a fundamental 
role in transport calculations for high-charge high-energy (HZE) ions and has been used to great 
effect by radiation investigators at the NASA, Langley Research Center. These earlier works 
have not however taken into account such important effects as straggling or of the energy 
downshift and dispersion which occur whenever a nuclear event takes place. In addition to the 
validation of physical processes, a theoretical model of the role of straggling is essential to 
understanding of the radiobiology of ion beams as required in evaluation of astronaut risks which 


must be minimized at least to within some regulated level [3]. Radiation therapy with ion beams 
also depends on the nature of the ion field near the end of the beam path in tissue [4] and is thus 
significantly affected by straggling. 

Achievements 

An analytic solution to the one-dimensional Boltzmann equation in the form of a Green’s 
function representing nuclear and atomic/molecular processes has been constructed. Based on the 
results of Wilson, Tweed, Tai and Tripathi [5] , a first order approximation to the ion beam 
range and energy straggling has been incorporated as a normal distribution for which the 
standard deviation is estimated from the fluctuation in energy loss events. Account has also been 
taken of the energy downshift and dispersion associated with nuclear events by the inclusion of 
appropriate terms in the nuclear transition cross sections (Tweed, Wilson and Tripathi [6]). A 
new computer code, based on this solution, is under development and some preliminary 
computations have already been performed. 

Some results from the new code have been compared with the earlier results of Wilson et 
al. ([7]), and were presented at a the DOE/NASA Radiation Investigators Workshop held in 
Washington, D.C., June 2001 [8], and at The World Space Congress , Houston, TX, October 
2002 [ 6 ], 

Parts of the code have also been used to develop a new model of a spacesuit liquid 
cooling and ventilation garment (LCVG). Past evaluations of spacesuit shielding properties 
assumed the basic fabric lay-up and LCVG could be homogenized as a single layer 
overestimating the protective properties over 60 percent of the fabric area. The new spacesuit 
model represents the inhomogeneous distributions of LCVG materials (mainly the water filled 
cooling tubes). An experimental test was performed using a 34-MeV proton beam and high- 
resolution detectors to compare with model-predicted transmission factors. The results of this 
study were reported at the above mentioned DOE/NASA Radiation Investigators Workshop[9]. 
and also at the 31 st International Conference on Environmental Systems, Orlando, Florida, July 
2001 [ 10 ]. 

Further space suit fabric studies were also reported at 32 st International Conference on 
Environmental Systems, San Antonio, TX, July 2002 [1 1] as were contributions towards the 
development of a new 3-D space transport code [11], 

Attachments 

The study funded by this grant has resulted in the publication of papers 5,6,10,1 1, and 12 
copies of which are attached below. 
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7 Abstract 

8 Ion beam transport theory allows testing of material transmission properties in the laboratory environment generated by particle 

9 accelerators. This is a necessary step in materials development and evaluation for space use. The approximations used in solving the 

10 Boltzmann transport equation for the space setting are often not sufficient for laboratory work and those issues are the main 

1 1 emphasis of the present work. In consequence, an analytic solution of the linear Boltzmann equation is pursued in the form of a 

12 Green’s function allowing flexibility in application to a broad range of boundary value problems. It has been established that simple 

13 solutions can be found for the high charge and energy (HZE) by ignoring nuclear energy downshifts and dispersion. Such solutions 

14 were found to be supported by experimental evidence with HZE ion beams when multiple scattering was added. Lacking from the 

1 5 prior solutions were range and energy straggling and energy downshift with dispersion associated with nuclear events. Recently, we 

16 have found global solutions including these effects providing a broader class of HZE ion solutions. 

17 © 2004 COSPAR. Published by Elsevier Ltd. All rights reserved. 

18 Keywords: Radiation risk; Boltzman transport equation; Ion beam transport; An improved Green’s function 


19 1. Introduction 

20 In space radiation transport, the energy lost through 

21 atomic collisions is treated as averaged processes over 

22 the many events which occur over even relatively small 

23 dimensions of most materials and is referred to as the 

24 continuous slowing down approximation. It is reasoned 

25 that the few percent energy fluctuation in energy loss has 

26 little meaning for ions of broad energy spectra and es- 

27 pecially in comparison to the many nuclear events for 

28 which uncertainties are still relatively large. In contrast, 

29 the laboratory testing of potential shielding materials 

30 uses nearly monoenergetic ion beams in which the 

31 interpretation of the interaction with shield materials 

32 requires a detailed description of the interaction process 

33 for comparison to detector responses (Schimmerling 

34 et al., 1986). The development of a Green’s function 
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approach to ion transport facilitates the modeling of 35 
laboratory radiation environments and allows for the 36 
direct testing of transport approximations of material 37 
transmission properties. For a number of years, this 38 
approach has played a fundamental role in transport 39 
calculations for high-charge high-energy (HZE) ions 40 
and has been used to great effect by radiation investi- 41 
gators at the NASA, Langley Research Center. These 42 
earlier works have not, however, taken into account 43 
such effects as straggling or of the energy downshift with 44 
dispersion which occur whenever a nuclear event takes 45 
place. In addition to the validation of physical processes, 46 
a theoretical model of the role of straggling is essential 47 
to understanding of the radiobiology of ion beams as 48 
required in evaluation of astronaut risks which must be 49 
minimized at least to within some regulated level (Shinn 50 
et al., 1999). The present development is in the context 51 
of an asymptotic expansion of the 3D Boltzmann 52 
equation, for which, the lowest order term is along the 53 
forward ray. Additional asymptotic terms are discussed 54 
in an earlier work (Wilson et al., 1991) and a related 55 
paper (Wilson et al., 2002a). 56 
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2. The Boltzmann equation 


The specification of the interior environment of a 
spacecraft and evaluation of the effects on the astronaut 
is at the heart of the space radiation protection problem. 
For some time investigators at The NASA Langley 
Research Center have been developing techniques to 
address this problem and an in-depth presentation of 
their work is given by Wilson et al. (1991) although 
considerable progress has been made since that publi- 
cation (Cucinotta et al., 1998). The relevant transport 
equation is the linear Boltzmann equation. The lowest 
order asymptotic term is the straightahead approxima- 
tion. With the target secondary fragments neglected, 
Wilson et al. (1991), this equation takes the following 
form: 

= £ [ c jk (E,E')<p k (z,E')dE' 

k>j J 

- Vj{E)<t>j{Z,E), z>Z’, (1) 

where <j>j(z,E) is the flux of ions of type j moving along 
the z-axis at energy E in units of MeV/amu and <Tj(E) 
and a jk (E,E') are the media macroscopic cross-sections. 
The <r jk (E,E ') represent all those processes by which 
type k particles moving in the z-direction with energy E 
produce a type j particle with energy E moving in the 
same direction. Note that there may be several reactions 
which produce a particular product, and the appropriate 
cross-sections for Eq. (1) are the inclusive ones. The 
total cross-section Cj(E) with the medium for each 
particle type of energy E may be expanded as 

aj(E) = af{E) + a%E) + a){E), (2) 

where the first term refers to collision with atomic 
electrons, the second term is for elastic nuclear scatter- 
ing, and the third term describes nuclear reactions. The 
corresponding differential cross-section is given as 


°AE,E') = YetfnWjW-E' + £„) 

n 

+ <rf(E)8 Jk 6(E-E) + ^ 
x exp 


\A2n)e jk 
{E + X jk -E'f 


2 4 


(3) 


where e n are the atomic/molecular excitation energy 
levels and where the collision energy downshift X jk and 
corresponding energy width e jk are approximated from 
the known momentum distributions observed in heavy 
ion reactions and represented by a gaussian model. 
Many atomic collisions (MO 6 ) occur in a centimeter of 
ordinary matter, whereas ~10 3 nuclear coulomb elastic 
collisions occur per centimeter, while nuclear reactions 
are separated by a fraction to many centimeters 
depending on energy and particle type. This ordering 


allows flexibility in expanding solutions to the Boltz- 
mann equation as a sequence of physical perturbative 
approximations. 

We require to solve Eq. (1) subject to a boundary 
condition of the type <j>j(zf,E) = Fj(E). In the case of a 
unit source at the boundary, Fj(E) takes the special form 

Fj{E) = 5 jk 5{E-E' ), (4) 

and the corresponding solution, which is called the 
Green’s function, is denoted by the symbol G jk (z,z ! , £,£*')• 
Once the Green’s function is known the solution for an 
arbitrary boundary condition Fj(E) is then given by 

e) = ^2 f E > e 'W) dE "- ( 5 ) 

In the case of an accelerator beam, the boundary condi- 
tion consists of a narrow gaussian function in energy and 
is incorporated by addition to the straggling width on 
leaving the boundary. In the case of space radiations, the 
boundary condition is represented as a broad function of 
energy arid direction for each ion type and is handled by 
ordinary numerical procedures. It should also be noted 
that Eq. (5) provides a basis for multiple layers of mate- 
rials by matching the solution at the boundary interface. 


3. Solution methods 


We rewrite Eq. (1) in operator notation by defining a 


vector array field function as 
4> = \4>j{z,E)}, 
the drift operator 


(6) 

D = [9,], 



(7) 

the interaction operator 



/ = S-<r = 

J <j jk (E,E')dE' 

- M*)], 

(8) 

with the understanding that / has three parts associated 
with atomic, elastic, and reactive processes as given in 
Eqs. (2) and (3). Eq. (1) is then rewritten as 

& 

li 

•e- 

Q 

= [/ at +r [ + / r ] • 


(9) 


and one must look for solutions. In what follows, we 
will recall the solution of the atomic interactions by 
Payne (1969) and implemented by Wilson et al. (2002b). 
Effectively, we look at 

£>•$=/*#, (10) 


which must then be coupled to the remaining terms in 
Eq. (9). For analysis, it will be advantageous to make 
the following separations: 


[D-I* — V +er r }-<P = 


j <T r jk (E,E')dE' 


<P = Z T <P. 
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142 3.1. Atomic processes 

143 The lowest order approximation to the Boltzmann 

144 equation is given in terms of the atomic collision pro- 

145 cesses as 

£>< P = / at (12) 

147 with the boundary condition 

* B = [ WM = [***{£-£)). (13) 

149 The solution, which incorporates energy straggling, 

150 takes the form 


<t>j(z,E) 


h 

V2ns' k (z - z 1 ) 


exp 


2 J’ 

(14) 


152 where 

(E k {z-z!)) = R;'[R k {E')-{z-l)}, (15) 

154 where R k (E) is the usual range-energy relation and 

155 /*(z-zO is the rms deviation for incident £-type parti- 

156 cles of energy E 1 after a distance of penetration z-z 7 

157 (Wilson et al., 2002b). 

158 3.2. Elastic scattering processes 

1 59 The addition of elastic scattering processes is given by 


D <P = [/“+/“] •#. (16) 

161 Since we have approximated the elastic scattering dis- 

162 tribution by 

<T'' k (E,E') = <T?(E')d Jt 5(E-E% : (17) 

164 we find that 

[/ cl M«[0], (18) 

166 and thus 

D -#«[/ at ] •$. (19) 


168 Elastic scattering does not appear in the first asymptotic 

169 term evaluated herein. The first correction will contain 

1 70 elastic scattering as a dominant term for the propagation 

171 of the surviving primary beam ions and in some bound- 

172 ary problems involving collimators elastic scattering will 

173 play a role for higher order terms. The elastic scattering 

174 propagator is a focus of current research and will couple 

175 with the present formalism. In the past, this coupling was 

176 in terms of acceptance functions and provided good 

177 agreement with neon ion beams (Shavers et al., 1993). 


[D - r +<r t }-4>= J ff' jk (E, E r ) d E' 


<P 


— 5 r . # 


( 20 ) 

In the present work, we approximate the fragment en- 
ergy distribution by 






(E + kjk — E) 

2* 


( 21 ) 


where X jk is the collision energy downshift (MeV/amu) 
and £ jk is the interaction energy w'idth (Mc\ /ainu^, A.j k is 
related to the momentum downshift (MeV/c) 


{11) 

via the equation 

i P(E)P* n3) 

** -AArn+E)’ { ’ 

where A k is the projectile mass (amu), Aj is the fragment 
mass (amu), E is the fragment energy (MeV/amu), m is 
the energy equivalent of a proton mass and 


p{E) = y/E 2 +2mE, (24) 

is the fragment momentum (MeV/amu/c). The interac- 
tion energy width is similarly related to the momentum 
width o F (MeV/c) through the equation 

p{E)or 


ejk ~Aj{m + E)’ 
where c r F is given as (Tripathi et al., 1994) 


<7p 


- l J 45 25 \( A^A* ~Aj) \ 

“\2 U* aT)\ 1 J' 


(25) 


(26) 


We start with the solution of the equation 

[Z>-/ at + a r ].G° = [ 0], (27) 

for a unit source at the boundary. Note that G° is di- 
agonal and takes the form 

jiK Pj{E) >/Stf t (z-*») 


x exp 


(E-wz-m 2 

2®i(z -z') 2 


(28) 


where the nuclear attenuation is described by the func- 
tion 




(29) 


178 3.3. Nuclear reactive processes 

179 Following the above analysis, we are left with 


and S k (E) is the change in E per unit path length per 
nucleon. Eq. (28) and the reactive integral operator are 
all that is required to develop the solution under the 
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straightahead approximation. The lateral spread of the 
beam is beyond the scope of the present development. 
So far all of the operators have had only diagonal ele- 
ments. Off-diagonal elements enter through the reactive 
regeneration terms <f jk which appear on the right side of 
Eq. (20). The challenge is to further develop the solution 
of Eq. (20) and this will be accomplished as follows. The 
integral form of Eq. (20) can be written as 

# = | D-I at + oT‘ • 4 >b + J [/>-/"' + (T r r‘ • 2 ■ <P dz, 

= <7° <p B + O ■ <7° • S* • <?, (30) 

where $ B is the appropriate boundary condition. Eq. 
(30) is a Volterra integral equation and is easily solved in 
a Neumann series as 

<? = [G <) + 2'G 0 'S r -G 0 + e'G 0 'S r '2‘G 0 -S r 
• Q ■ G° + *■•]* #b 

= [G° + G x +G 1 -{ ] • # B , (31) 

with the elements of the leading term given as Eq. (28). 
The above formalism lends the following interpretation 
of the solution. The operator G° propagates the particles 
with attenuation processes. The first term <j° • # B 
propagates the ions at the boundary to the interior. 
E x • G° * # B is the production density of first generation 
secondaries at depth z\. These are propagated to the 
interior by • E x • G° ■ Lastly, G 1 • = Q 

G° • S* * G° • # B represents the sum of all the first gen- 
eration secondaries being propagated from the interval 
[z!,z] and so on. We have already identified the propa- 
gator G°. We now need to identify the remaining terms 
in the Neumann series and we begin noting that these 
are related via the recurrence formula 

G fl+1 =[0.G°.S r )-G", n>0. (32) 


3.4. First collision term 

The second term in Eq. (31) is the first collision term 
G) k (z,J,E,E?) = IQ- G° ■ S' • 

= jf J G°(z,z,.£.£ 2 ) 

* {/ 

xG°*(z,z\£,£')d£i}d£2dz,. (33) 

The physical interpretation is that E x * G° is the volume 
source of ions from collisions at z\ of a unit ion source at 
z / of energy E f The ions present at z with energy E are the 
result of propagation from the all the ions through out 
the volume. The first task is to evaluate the volume 
source term 


[S’ ■ G°] t (zi,z' ,E 2 ,E') 


f •}.«> „„ 

(E 2 + A,t - Ei) 2 ] 

1 v!*!* “ p 

<N 

t 

i 


, ft(g) 1 

Pk(Ei) \/2ns' k (z - z 1 ) 


exp 


(£,-<^(z-z'))) 2 
2 s' k {z-zff 


d£,. 

(34) 


Note that a sharp maximum occurs at £i = (E' k (z - z')), 
E 2 = E\- A;* and the cross-sections and attenuation 
functions are siowiy varying functions of energy so that 
Eq. (34) can be accurately approximated as 


[S r G%(z,,z',£ 2 ,£') 

- iw&g W*-™ 


^ 27 l[si(z-2') 2 +4; 


xexn / (£2 + / M ~ (^( Z - Z 7 ))) 2 1 

P l 2 [s' k (z - z 1 ) 2 + e} k ] j 

The next step is to construct the term 


(35) 


[<7° • S' • Cp)j k {z,z\,d ,E,E?) 

pm 


-h 


exp 


(E)V2ns’j(z — Z]) 

Pk(E') °jk(( E 'k( z i -*))) 


f (£-(£»(z-z,))n 

I 2Sy (z — Zl ) 2 / 




xex 

' 2[s> k (z-z>) 2 + E %\ 


jd£ 2 , 


(36) 


where (E"(z - z x )) = RJ 1 [Rj(E 2 ) - (z - z,)] and 
s"(z - zi) is the corresponding spread. The integral has a 
sharp maximum at E 2 = {E[{z\ - z!))k jk = {E 2 ), E — 
( E"(z-z\ )) = ( Ej(z )), where the cross-sections, attenua- 
tion functions, and straggling widths are evaluated. We 
expand ( E"(z-z\ )) about the maximal value of E 2 to 
obtain 


z(£)'(z -z,)) « (Ej(z)) + r jk [E 2 - ((E' k (z, - /)) - A,,)], 

(37) 

where 

rjk = [6 £2 (£" (z - z, ))] <£2> = - (38) 

Substituting Eq. (37) into the integral (36), making 
the change of variables x = r Jk \E 2 - {(E' k {z\ - z!)) — A,*)] 
and integrating with respect to x results in 


250 

251 

252 

253 


255 


257 

258 

259 

260 
261 
262 
263 


265 


267 

268 
269 



JASR 6703 
TYPED/ 3/4/04 


ARTICLE IN PRESS 


J. Tweed et al / Advances in Space Research xxx (2004) xxx-xxx 


No. of pages: f 
DTD 4.3.1 / SPS-N 


[G° • S' • G |0 ]^(z,z 1 ,z' l £,£') 

Pj{E)P k {{E' k {z x -*))) 

» z£)l rrn { \ E zlM) \ 
^0 l J’ 

271 where 


(39) 


/y(zi) = *,•' {^[(£i(z. - z')) - *jk\ - (z - z,)} (40) 

273 and 

•V*(*i) = ^ ^ [ii(zi - ^) 2 + 4] + <(z - Z| ) 2 - (41) 

275 Lastly, we need to evaluate the integral 

G) k {z,z!,E,E) = {Q(?-S'- G°]y*(z, zi,^,E,E) 

= J\(?-3'-(?} jk (z,z l ,z',E,E')dz l . 

(42) 


277 For a given set of parameters z, z\ there is a value 

278 z m of z\ at which the integrand of (42) achieves a max- 

279 imum and at which slowly varying factors entering the 

280 integrand may be computed. Thus 


282 

283 


G) k (z ) z',E,E f ) 


P J ((E' k (z m -z / ))->.j k )P k (E) 
Pj(E)P k ((E' k (z m -z '))) 

> , ^«(£;(z„-z'))) 

\l 2 nsj k {z m ) 


X 



|E~/;(*i)] 2 

2Sjk( z m) 


jdz,. (43) 


The point z m at which the integrand of (43) achieves its 
maximum is given by the equation 


fj{z m )=E, 


(44) 


285 and is easily obtained by the routine root finding tech- 

286 niques. It is not difficult to show that 


P.(z 1 = Ji &[(£&,-*'))] 
JjKm) 1 Sj[(E' k {z m -zf)) - kj k \ 



(45) 


288 Therefore, in (43) we may use the substitution 

289 x = [£ -/y( Z |)]/[\/5^i(z m )] and then integrate to get 


G) k {z,z',E,E) 


PjKE^-^-XME} 
Pj(E)P k ((E' k (z m — z 1 ))) 

a%m(z m -/)))) 

2 fj(z m ) 


< erf 


— erf 

' E — fj(z ) " 

t 

_VlSj k {z m )_ 


_V2Sjk{Zm). 


} 


( 46 ) 


3.5. Second collision term 292 

The third term in Eq. (31) is the second collision term 293 
G\{zJ,E,E' ) = [Q- G° • S' • G%{z^,E,E) 

= J‘ J tfj(z,z u E,E 2 l) 

* { 

l p>j J 

x Gj rf .(zi,z',E 1 ,£')d£i|d£ 2 dz|. (47) 

On making approximations similar to those used in the 295 
previous section, Eq. (47) is reduced to the form 296 


z.om 



Fig. 1. Primary ion flux at various depths for Ne(20,10) incident on 
aluminum at 600 MeV/amu. 



291 where s,*(zi) is given by (41). 


Fig. 2. First generation 0(16,8) fragment flux at various depths com- 
pared with previous results (broken curve). 
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G) k {z,£,E,E') Pj[E]s\e]^ Ej + 

x G 'p* ( z i > z '» E i + E')iz x , (48) 

298 where 

Ej=R-'{Rj{E)+z-zi\, (49) 

300 and is then evaluated by numerical quadrature. 


301 4. Results 

302 Shown in this section are some results for the gO 16 

303 fragments which are produced when a beam of ioNe 20 



Fig. 3. Second generation 0(16,8) fragment flux at various depths 
compared with previous results (broken curve). 


ions strikes an aluminum target at 600 MeV/amu. The 304 
results presented are similar to those obtained for other 305 
ions. 306 

Fig. 1 shows the flux of the primary beam at various 307 
depths and exhibits the effects of energy straggling. In 308 
contrast to earlier works in which the primary beam 309 
appears as a propagating delta function, we see here that 310 
the primary beam attenuates and widens with depth. 311 
Note that the greatest depths in Fig. 1 are beyond the 312 
85% range where straggling propagators of the past 313 
have failed. 314 

Figs. 2 and 3 show the flux of the first and second 315 
generation of gO 16 ions produced, respectively, in com- 316 
parison with the corresponding flux (broken curves) 317 
obtained from an earlier approximate code using a non- 318 



Fig. 5. Second generation 0(16,8) fragment flux at various depths 
compared with the same flux for the case in which the collision energy 
downshift is zero (broken curve). 




Fig. 4. First generation 0(16,8) fragment flux at various depths com- Fig. 6. First generation 0(16,8) fragment flux at various depths com- 
pared with the same flux for the case in which the collision energy pared with the same flux for the case in which the interaction energy 

downshift is zero (broken curve). width is zero (broken curve). 
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319 perturbative expansion of the solution (Wilson et al., 

320 1991). The non-perturbative approximation lacks spec- 

321 tral details and assumes a broad near uniform distri- 

322 bution over the allowed energy domain (Wilson et al., 

323 1991). 

324 The effect of the collision energy downshift Xjk is ex- 

325 hibited Figs. 4 and 5, where the flux of the first and 

326 second generation of gO 5 * * * * * * * * * * 16 ions is compared with the 

327 corresponding results for the case in which the down- 

328 shift is zero. For the ions shown the shift is not great, 

329 contributing only a few MeV/nucleon. The downshift of 

330 more massive projectiles is somewhat larger. 

331 Figs. 6 and 7 exhibit the effect of the interaction en- 

332 ergy width £,* by comparing the flux of the first and 


Fig. 9. Second generation 0(16,8) fragment flux at various depths 
compared with the same flux for the case in which there is no energy 
straggling (broken curve). 


second generation of gO 16 ions with the corresponding 333 
results for the case in which the interaction energy width 334 
is zero. Significant widening occurs in the first genera- 335 
tion of secondaries but little effect is seen in the second 336 
and presumably higher generations. 337 

The effect of energy straggling on the the first and 338 
second generation of gO 16 ions is exhibited Figs. 8 and 9, 339 

where the flux of of these ions is compared with the 340 
corresponding results for the case in which no energy 341 
straggling is present. In contrast with the results for the 342 
primary beam where straggling makes a significant 343 
contribution the effect on the first and second generation 344 
of secondaries appears to be relatively small. 345 




E, MeV/nudeon 



E, MeV/nudeon 


Fig. 8. First generation 0(16,8) fragment flux at various depths com- 
pared with the same flux for the case in which there is no energy 
straggling (broken curve). 


5. Concluding remarks 346 

The present formalism provides means of easy vali- 347 

dation of material transmission properties in conven- 348 

tional laboratory setups at least to the third 349 

perturbation term. Higher order terms can be easily 350 

added using non-perturbation theory and assuming the 351 

third term spectral distribution. The next step is the 352 

simulation of the detector responses so that “raw” ex- 353 

perimental data can be used to validate model predic- 354 

tions thereby simply fying the validation process. Early 355 

versions of the Green’s function code, when coupled 356 

with multiple elastic scattering in terms of acceptance 357 
functions (Shavers et al., 1993), showed great promise in 358 
describing HZE ion transport. In those studies, the 359 
straggling, energy downshift, and dispersion were ne- 360 
glected. The present formalism corrects those last re- 361 
maining deficiencies. The recognition of the present 362 
formalism as the lowest order asymptotic term provides 363 
a systematic approach to more realistically treat a host 364 
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of ion beam related problems. The next step will be to 
couple the multiple scattering propagator to the for- 
malism and adding transverse momentum components 
to the first interaction term. 


6. Uncited references 

Schimmerling et al. (1999); Tschalar and Maccabee 
(1968). 
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